Efimov Physics and the Three-Body Parameter within a Two-Channel Framework 
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We calculate shallow three-body bound states in the universal regime, defined by Efimov, with 
inclusion of both scattering length and efi'ective range parameters. We find corrections to the 
universal scaling laws for large binding energies. For narrow resonances we find a distinct non- 
monotonic behavior of the threshold at which the lowest Efimov trimer merges with the three- 
body continuum. The origin of the three-body parameter is related to the two-body atom-atom 
interactions in a physically clear model. Our results demonstrate that experimental information 
from narrow Feshbach resonances and/or mixed systems are of vital importance to pin down the 
relation of two- and three-body physics in atomic systems. 
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Introduction. The counterintuitive behavior of three- 
body systems at the threshold for two-body binding is 
highlighted by the Efimov effect where an infinitude of 
geometrically scaling states appears [U. While unsuc- 
cesfully sought for in nuclear physics |2'|, the effect has 
been confirmed and explored in ultracold atomic gases 
Q . From the theoretical point of view these systems have 
been described accurately by universal theories that only 
take the lowest order scattering dynamics into account 
through the two-body scattering length, a d, @ . How- 
ever, the overall scale of the spectrum cannot be obtained 
in the universal theory and the so-called three-body pa- 
rameter, A, is needed to complete the formalism. 

In ultracold atomic gases, trimer physics can be stud- 
ied using inter-atomic Feshbach resonances @ that pro- 
vide tunability of a over many orders of magnitude. The 
parameters of these resonances will in general depend 
on the microscopic details of a particular atomic system. 
Still, in a surprising development, the Grimm group have 
reported strong indications that the three-body parame- 
ter is determined by the van der Waals length, rvdw 0- 
This implies that atomic trimer physics in the weakly 
bound limit only depend on two-body parameters. The 
experimental findings have generated a flurry of recent 
theoretical interest [sl-lllj. It has been suggested that 
the presence of many deep bound states in the two-body 
potential typical of alkali atom systems will suppress the 
dependence on short-range physics due to a large inner 
repulsive barrier [§, [T^f . 

The most natural way to eliminate any need for a 
three-body parameter is to use a finite-range two-body 
potential that exhibits the known features of the inter- 
atomic interaction. Here we pursue this natural and 
straightforward approach to the issue of the origin of A 
in three-body physics. As three-body calculations with 
finite-range potentials are typically very cumbersome, it 
is highly desirable to have simple models that mimic the 
finite-range physics. Therefore we consider here a model 
with zero-range interactions that includes finite-range 
correction effects. Consequently, we solve the three- 
body Schrodinger equation using one- and two-channel 
models for the description of the Feshbach resonances 
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FIG. 1: (Color online) The threshold scattering length a'~' 
at which the lowest universal Efimov trimer merges with the 
three-atom continuum plotted against the strength, Sros, of 
the Feshbach resonance. The right-hand side corresponds to 
broad resonances. The different curves of like linetype show 
result with different three-body parameters, pc, in units of 
the atomic van der Waals length, rvdw- The values of pc 
decrease from top to bottom (precise values are given in the 



text). Experimental values are from ^^^Cg ^ 
[13, and ^'^Rb [Hi. 



that are used to control interactions in experiments [6j. 
The three-body parameter is considered parametrically 
in both models, and this allows us to study the thresh- 
old behavior of Efimov trimers as function of A and the 
width of the Fesbhach resonance which is related to the 
effective range, rg Q. Here rvdw and r^. are intimately 
connected through semiclassical calculations [IMl. It 
is conceivable that this explains the observed values of 
A as Te or Tvdw provides a background scale that de- 
termines the overall scale. However, given the delicate 
non-classical nature of the universal trimer states it is 
far from obvious if and how this can work out. 

The experimental findings indicate that for several dif- 
ferent atomic species the ratio of threshold scattering 
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length a^~^ for creation of Efimov trimers out of the 
three-body continuum is of order a^~'' /r^dw = 8-5 ~ 9.5 
0. Here we show that this result can be obtained in a 
single-channel zero-range model for a specific choice of A, 
which can be related to the underlying two-body atom- 
atom physics in a natural way. Within a two-channel 
model we find an intriguing non-trivial behavior of a^") 
for narrow resonances, irrespective of the three-body pa- 
rameter, see Fig. [TJ We also study the behavior of a^") 
as a function of the number of bound states allowed by 
the two-body atomic potential for both one- and two- 
channel models. Generally, we find that the inclusion of 
effective range decreases |a(~^|. Our results predict that 
narrow resonance systems are important for obtaining a 
full picture of the relation between two- and three-body 
parameters for universal bound state physics. 

Method. We consider a system of three identical 
bosonic particles using hyperspherical coordinates fl^ 
defined from the Cartesian coordinates ri,rj,rk through 
Xi = (rj-rfe)/\/2 and yi = |(ri - (r^ +rfc)/2) as hyper- 
radius, p = \/x^ +yf, and hyperangle, ai = tan"^ 
{i,j,k} are cyclic permutations of {1,2,3} and p is in- 
dependent of this choice. We apply the hyperspheri- 
cal adiabatic approach with wave function il) — 
P~^''^T.nfn{p)^n{p, ^) , whcrc = {a,, X, / \x,ly, / \y,W 
is a set of angular coordinates. We keep only the lowest 
adiabatic potential corresponding to n = 0, and index n 
is henceforth suppressed. For the description of Efimov 
trimer states, this approximation has proven extremely 
accurate 3. The radial equation is 

where \{p) is the eigenvalue to the hyperangular equation 

(^A2 + ^l/)cf>(p,17) = A(p)ci>(p,17), (2) 

in which is the generalized angular momentum oper- 
ator, V is the two-particle interaction potentials, and m 
is the atomic mass. In Eq. ([T]), the non-adiabatic correc- 
tions are omitted as they are found to be negligible. 

We use zero-range single-channel and a two-channel in- 
teraction models [2y]. The former is accurate for broad 
Feshbach resonances (small r^) and the latter for narrow 
resonances (large r-g). Both models contain the back- 
ground scattering length in the open channel, flopcn- Our 
model must predict three-body properties solely from 
two-body potentials, so we must insist that flopon is a 
quantity most naturally associated to the two-body rvdw- 
In a semiclassical approach, the scattering length of such 
potentials is a = 47r(r[l/4])~^rvdw ~ 0.956rvdw [Hi- 
Then Oopcn should be identified with Tvdw to within a 
few percent. The resonance strength is related to re [20| . 
We characterize the strength by Sros = ^Vdw/I-Roli where 
i?o is the effective range when \a\ = oo [a]. 

Eq. ID) is solved numerically for /(p) with the con- 
dition /(pl) = for large pL- At short distance the 




FIG. 2: (Color online) Trimer binding energy Et vs. scatter- 
ing length a for the one- and two-channel models with effective 
range re = — 5rvdw. is the threshold scattering length 

for existence of the lowest trimer. Finely dotted lines indicate 
atom-dimer threshold for positive a. Successive states have a 
scaling of 515 1] and for visibility the axes are scaled by the 
power 1/8. 



zero-range models used here require a cut-off, pc, with 
f{pc) = [ij. Initially, we consider pc a parameter and 
study the three-body spectrum for different pc and dif- 
ferent Sros- Pc is the coordinate-space equivalent of A. 
Below we relate the pc to the two-body atomic potential. 

Results. In Fig. [5] we show trimer energies, Et^ as 
function of a for both one- and two-channel models. We 
find that the two-channel energies are generally lowest 
when a < and for a > the trimers cross before 
merging with the atom-dimer continuum. As indicated 
in Fig. [21 a'-^'' is the threshold scattering length for ap- 
pearance of the lowest Efimov trimer on the negative a 
side. The two-channel model is here seen to move this 
threshold to the left, i.e. to smaller negative a^~\ These 
thresholds and the energy for \a\ = oo are connected 
by universal relations, where lower energy on resonance 
translates to smaller |a(~)| at the threshold [1]. 

A systematic study of the influence of both pc and r^. 
(or equivalently Sros) is shown in Fig. [1] which is one of 
our main results. The values of Pc/rvdw from top to 
bottom in Fig. [T] are 1.20, 0.82, 0.66, 0.58, 0.51, 0.47, 
0.42, 0.40, and 0.38. They correspond to n = to 9 
in Eq. ^ below. Both models agree for s^cs ^ 1 and 
we plot the two-channel model results only in the region 
where it deviates. We find that to reproduce the exper- 
imental data for s^es ^ 1, a cut-off of Pc/fvdw — 0.58 is 
required. However, for small s^es, the same cut-off does 
not reproduce the known data point coming from ''Li 
(other measurements have smaller a^~^ [SI, H^], which is 
closer to our predictions). While we do find an increase 
toward the "^^K data point at small s^es, it cannot be 
accomodated for the same pc. 

In general, we find that, irrespective of pc, the inclu- 
sion of effective range brings a non-monotonic behavior 
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to a^~\ and it tends to push the value of a^~^ down 
for smah Sj-cs or large |re|. The reason can be seen m 
Fig. |6] where the adiabatic potentials for different effec- 
tive ranges, r^/rvdw = 0, — 0.1, — 1, — 5, are shown. At 
rv = the bound states generally reside at large p/r^d\N 
[J|. The effective range causes an additional repulsive 
barrier initially leading to less bound energies when |re| 
is small, since the bound state sitting at large p will feel 
the barrier at first. When the effective range increases 
the attractive region will ultimately lead to an increase 
in bound state energy as the wave function starts to oc- 
cupy the pocket at smaller p. 

This is in sharp contrast to the study in Ref . jllj which 
finds the opposite behavior. As argued above, it is phys- 
ically reasonable since Et is lower at finite re < 0. We 
notice that our two-channel model has rg < which is 
consistent with the usual theory of Feshbach resonances 
@. Ref. [ll| appears to accomodate also re > and this 
may resolve the discrepancy. 
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FIG. 3: (Color online) The lowest adiabatic potentials multi- 
plied by as functions of hyperradius, p, for the zero range 
one-channel model and three different effective ranges for the 
two-channel model. 

Connection to two-body potential. The zero-range 
model studied above does not carry information about 
the van der Waals length as it stands. However, the 
three-body parameter or cut-off, pc, has a physical mean- 
ing as it provides a hard-core repulsion in hyperspherical 
three-body coordinates. To connect the formalism to ex- 
perimental data, it is therefore necessary to find a rela- 
tion between the two-body atomic physics and pc- Below 
we consider different two-body interaction models and 
derive analytical formulas for the behavior of a'^' as the 
two-body interaction parameters are varied. 

In Fig. 13] we plot Morse, Lennard- Jones, and a van der 
Waals potential with a hard-core at Tc- We first focus on 
the van der Waals plus hard-core model. 
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FIG. 4; (Color online) Typical two-body neutral atom-atom 
potentials, including the hardcore van der Waals potential of 
Eq. the Lennard- Jones potential, VLj{r) — — and 

the Morse potential, Viviir) = 4eo(e~^"<''~''«' - e""'''"''"'). 
Here Ce, Cio, and a are parameters. They have all been 
fitted to have the same overall energy, eo, at the radius of the 
potential mimiiiiiiii, /"mill- 



where Tvdw = (^7^)^^^- In order to relate the behav- 
ior of this potential to the physics of a Feshbach reso- 
nance, we use the formula |16l |l8| a = a(l — tan <I>) with 

•I* = —f¥^ — ^ where a diverges when $ = {n + ^)tt 
for integer n, which counts the number of s-wave bound 
states accomodated by the potential. Thus 



2_ / rvdw\ _ 7 
tt\ Tr) 8 



(4) 



rounded to highest nearby integer. 

In the zero-range single-channel model pc and a^") are 
related, since pc is the only length scale available. How- 
ever, it is a non-trivial matter to determine this relation 
within the model. Numerically, we find the linear rela- 
tion a(~) = —5 Pc with S = 31.756. To cast this relation 
into a form that depends only on two-body physics, we 
observe that the two-body hard-core at Tc is also respon- 
sible for a three-body hard-core cut-off, more precisely 
Pc — V2rc- This condition ensures that when each two- 
body subsystem has radius r > Tc, the third particle will 
also be outside with respect to the others (see Ref. [23j 
and the supplementary note|X]). We note that this rela- 
tion is not the same as the pc used in Ref. 9]. Using 
a^^-* = —5 Pc we obtain 



7.(-) 



25 



(5) 



r > Tc 



(3) 



This semi-analytical expression for the threshold in terms 
of the number of bound states is one of our main results. 

The relation in Eq. ([S]) is plotted in Fig. [5] along with 
the experimental data and the numerical results obtained 



4 



25 
20 
15 
10 

5 



» 39j^ Mean value 

One-channel 
Two-channel, log(srcs) = 
Two-channel, log(srcs) = —0.7 




10 



15 



20 



25 



FIG. 5: (Color online) Semianalytic results for a'~' plotted 
against number of bound states in the two-body van der Waals 
plus hard-core potential. The horizontal position of the ex- 
perimental data is arbitrary. The grey band indicates the 15% 
margin around the mean value of ~ 9.8. 



from the two-channel model for different value of s^cs- 
The one-channel model is consistent with data for n ^ 
10 — 20 and fits the universal ratio of Ref. [o'l for n = 
13. This is also consistent with the findings of Ref. [lo| . 
although the data only goes to n = 10. Actually, the 
nT^f^ behavior that we find here seems to also appear in 
Ref. [lo| . where an extension to higher n could confirm 
this prediction. 
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Re [A] 


2.67 


3.92 


4.18 


4.65 


De [eV] 


1.06 


0.52 


0.49 


0.45 


riLj 


41 


99 


152 


201 


riM 


28 


67 


103 


137 



TABLE 1: Potential parameters, (bond length) and 
(dissociation energy) [S^l, and estimated number of bound 
states for Lennard- Jones, ulj, and Morse, um, potentials. 

Even more interestingly, our results for small Sres h> 
dicate that \a^~^\ drops faster with n than for Sres ^ 1- 
This is seen in the experimental data on ^Li which is 
slightly below the ^^Rb and ^■^'^Cs points, but our model 
seems to overestimate this trend. Clearly, more results 
on narrow resonance systems are required to address the 
question of effective range corrections. We expect a lower 
ja^^^l value than for broad resonances. 

Two-body potential models Above we employed a van 
der Waals plus hard-core potential. However, as seen in 
Fig. m more realistic Lennard- Jones or Morse potentials 
have a smooth behavior of the inner barrier. This implies 
only minor quantitative corrections that nevertheless de- 
serve to be addressed along with the number of bound 
states expected in a real alkali dimer. 



The number of s-wave bound states in the Lennard- 
Jones (LJ) and Morse (M) potentials can be estimated 
analytically and yields [231 '^iJ = O.SGly^— | and um = 

0.245v^ — 5, where /3 — '"2"^'^° with Tmin the radius 
at which the potential takes its minimal value, eg. For 
comparison, the expression in Eq. ^ can be written n — 
0.225v^ — |, with rmin Tc in /3. The similarity of 
these expressions makes it clear that the behavior seen 
in Fig. [5] is generic and does not depend on our choice of 
two-body potential. The difference in constant provides 
only a minor quantitative change in the numbers. 

An important question, however, remains about the 
number of bound states in a real alkali dimer system. Us- 
ing the molecular dissociation energy, eg = D^., and the 
bond length, r,nin = Re, of Ref. [IJl, we list the numbers 
for Li, K, Rb, and Cs dimers in Tab. H] The estimated 
number of bound states is outside the axis in Fig . [5] and 
also much beyond the results shown in ReL [lOj. The 
agreement with theory at a rather limited number of 
bound states (n '-^ 10 — 20) is then quite surprising. 

A number of important observations can be made. 
First, the decrease of la^"-*! with n is weak, and a shift of 
the length scale in Fig. [S] would therefore place the one- 
channel model within the experimental range for larger 
n and it would stay within the 15% deviation from the 
mean for a larger interval (since the slope at larger n 
decreases even faster). Second, the experimental data 
might indicate that only a certain number of bound states 
play an active role. Equivalently, even if the two-body 
potential is very deep, only the upper part of the two- 
body potential and the bound states closest to threshold 
set the scale of the three-body problem. This appears 
to be very reasonable since we are considering universal 
Efimov trimers here and not strongly bound three-body 
states. Third, the case of small Sj-cs has |a*^~''| oc 
with r > 1/2 as seen in Fig. [SJ This implies that narrow 
resonance systems should be even more insensitive to n 
beyond a certain lower limit. 

We can give a quantitative argument for the lack of 
sensitivity to the many deeper bound states in a van 
der Waals potential. The number of bound states, n{E) 
as a function of energy, E, counted from the E = 
threshold and down within the WKB approximation is 
n{E) = no{\E\/Eydwy^^, where no is the total number 
(given in Tab. U for different potentials) and i?vdw = 
ti^/nir^^yj- > is the depth of the potential [2^. For the 
Sros > 1 cases (^^Rb and ^^^Cs), n{E)/no ~ 0.10 - 0.20 
which implies |E|/i?vdw ^0.001-0.01. Numerically we 
find a three-body energy on resonance E^ = 0.006i?vdw 
(using pc = 0.58rvdw)- However, universality relates 
E3 — ti^K^/m and ~ —1.51 The energy 

scale at the continuum threshold is given by a^"^ through 
\E\ ~ ft^/m(a(^))^ — 0.003i?vdw, in agreement with the 
interval above. In the case of ''Li, E3 is similar but this is 
compensated by a smaller uq so this case can also be ex- 
plained. For heavier alkali atoms with narrow resonance, 
our two-channel results predict a smaller a^~^ /?'vdw than 
9.8, which is a good experimental test of our theory. 
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a) b) 

Triangular configuration Linear configuration 



FIG. 6: (Color online) Schematic of the a) triangular and 
b) linear configurations for an equal mass three-body system. 
ri2 and ri2,3 indicate the two Jacobi relative vectors. 



Conclusions. We derive an analytical formula that 
connects a^~^ and the number of s-wave bound states, n, 
in the two-body potential which agrees with recent ex- 
periments 71 ~ 10 — 20. While alkali atoms typically have 
larger n, we argue quantitatively that only a subset de- 
termine the properties of Efimov three-body states. We 
find non-monotonic behavior of a^~^ with increasing ef- 
fective range which demonstrates that systems with nar- 
row Feshbach resonances will provide both qualitative 
and quantitative understanding of universality and the 
origin of the three-body parameter. This means that 
different atomic species which typically have narrow res- 
onances and a denser spectrum would be very helpful. 
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Appendix A: Determination of the hyperradial 
boundary condition 

The atom-atom two-body potential has a steep repul- 
sive inner core which we will assume for the moment can 
be representated by a hard wall, i.e. an infinite potential 
for radii r < rc- In this case the boundary condition is 
simply that the two-body wave function must be zero at 
Vc- This must then be translated in the three-body prob- 
lem where it implies that the total wave function must be 
zero whenever any of the relative distances between two 
out of the three atoms is less than or equal to r^. Any 
penetration of the wave function into the wall would cost 
an infinite amount of energy and is thus forbidden. 

A neat and elegant way to obtain a condition on p is 
the following. As is easily shown the hyperradius fulfills 



n2 

rk) 



1 



1 



12 



12,3' 



(Al) 



where ri2 = ri - r2 and ri2,3 = ~ (ri -I- r2)/2 are 
respectively the relative vector from particles 2 to 1 and 
the relative vector from the center of mass of 1 and 2 to 
particle 3 (see FiglH]). 

Now consider the triangular configuration shown in a) 
of Fig. ini From the formula in Eq. (|Aip . it is clear that 
when ri2 — rc this setup also has p — Vc- However, the 
linear configuration of b) in Fig. |6] with p — rc will then 
have two atoms that are within the hard-core radius and 
yield an infinite contribution to the potential energy (this 
happens f.x. when \x\ = rc and y = putting atom 3 
midway between 1 and 2). 

Consider instead the linear configuration with ri2 — rc 
and impose the requirement r23 > rc, where r23 = r2 — 
r3. Since we have ri2,3 = r23 + ri2/2 we get 



2 \ 2 

P >rc 



2rt 



(A2) 



i<k 



We thus see that the condition p > y/2rc ensures that 
both triangular and linear configurations are outside the 
hard-core regions. Since these configurations are ex- 
tremal, the condition implies that no regions with infi- 
nite potential are sampled by the hyperradial three-body 
wave function. 

The rigourous formal argument for the validity of the 
relation pc — \firc using the hyperspherical approach 
can be found in Ref. [1^, where the relation is derived 
using a square well potential. The asymptotic region is 
precisely p > \f2rc as found above. Here we assumed a 
hard-core potential for simplicity which gives the factor 
of a/2. For a real atom-atom potential, the hard-core is 
slightly softer (typically of the which may lead to 

a minor change in the factor. 

The hyperspherical approach and the use of the low- 
est hyperradial potential to describe Efimov states is a 
well-known and accurate procedure. Additionally, as we 
have shown above, a hard-core repulsion in the two-body 
potential naturally generates a hyperspherical cut-off ra- 
dius, Pc = \/^rc. This constitutes a clear and simple 
explanation of the repulsive core found in the numeri- 
cal calculations in Ref. [13], and our results demonstrate 
that a simple model can accurately describe those nu- 
merically involved findings. We thus see how the value 
\/2rc for the three-body cut-off parameter naturally de- 
rives from the hard-core behavior of the corresponding 
two-body potential. 



[1] V. Efimov, Yad. Fiz 12, 1080 (1970); Sov. J. Nucl. Phys. [2] A. S. Jensen, K. Riisager, D. V. Fedorov, and E. Garrido, 
12, 589 (1971). Rev. Mod. Phys. 76, 215 (2004). 



6 



[3] F. Ferlaino and R. Grimm, Physics 3, 9 (2010). 

[4] E. Nielsen, D. V. Fedorov, A. S. Jensen, and E. Garrido, 

Phys. Rep. 347, 373 (2001). 
[5] E. Braaten and H. W. Hammer, Phys. Rep. 428, 259 

(2006). 

[6] C. Chin, R. Grimm, P. S. Julienne, and E. Tiesinga, Rev. 

Mod. Phys. 82, 1225 (2010). 
[7] M. Berninger et al, Phys. Rev. Lett. 107, 120401 (2011). 
[8] P. Naidon, E. Hiyama, and M. Ueda, Phys. Rev. A 86, 

012502 {2QV2y. 

[9] C. Chin, arXiv: 1111.1484 v2. 
[10] J. Wang, J. P. D'Incao, B. D. Esry, and C. H. Greene, 

Phys. Rev. Lett. 108, 263001 (2012). 
[11] R. Schmidt, S. P. Nath, and W. Zwerger, 

larXiv:1201.4310k 2 
[12] T. Kraemer et al, Nature 440, 315 (2006). 
[13] S. E. Pollack, D. Dries, and R. G. Hulet, Science 326, 
1683 (2009). 

[14] M. Zaccanti et al, Nature Phys. 5, 586 (2009). 

[15] R. J. Wild, P. Makotyn, J. M. Pino, E. A. Cornell, and 

D. S. Jin, Phys. Rev. Lett. 108, 145305 (2012). 
[16] G. F. Gribakin and V. V. Flambaum, Phys. Rev. A 48, 

546 (1993) 

[17] B. Gao, Phys. Rev A 58, 4222 (1998). 



[18] V. V. Flambaum, G. F. Gribakin, and C. Harabati, Phys. 
Rev. A 59, 1998 (1999) 

[19] Y. Suzuki and K. Varga, Stochastic Variational Approach 
to Quantum- Mechanical Few- Body Problems, (Springer- 
Verlag, Berlin Heidelberg, 1998). 

[20] P. K. S0rensen, D. V. Fedorov, and A. S Jensen, Few- 
Body Systems, 0177-7963 (2012) 

[21] N. Gross, Z. Shotan, S. Kokkelmans, and L. Khaykovich, 
Phys. Rev. Lett. 103, 163202 (2009). 

[22] N. Gross, Z. Shotan, S. Kokkelmans, and L. Khaykovich, 
Phys. Rev. Lett. 105, 103203 (2010). 

[23] A. S. Jensen, E. Garrido, and D. V. Fedorov, Few-Body 
Syst. 22, 193 (1997). 

[24] Igel-Mann et al, J. Chem. Phys., 84 9 (1986) 

[25] G. D. Mahan, M. Lapp, Phys. Rev., 179 1 (1969) 

[26] The number of bound states, n{E), can be found 
within the WKB formalism by considering the integral 
JJ° \J E — V(r)dr — -Kh{n{E) -\- 1/4) which is a measure 
of the number of nodes in the wave function with energy 
E < 0. Using E — one obtains the total number, no. 
The van der Waals potential does not yield analytical re- 
sults, but we numerically find that n{E) oc |-E|^/^, which 
when normalized yields the relation in the text. 



